Analytical properties of horizontal visibility graphs in the Feigenbaum scenario 
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Time series are proficiently converted into graphs via the horizontal visibility (HV) algorithm, 
which prompts interest in its capability for capturing the nature of different classes of series in 

■ a network context. We have recently shown [1| that dynamical systems can be studied from a 
novel perspective via the use of this method. Specifically, the period-doubling and band-splitting 
attractor cascades that characterize unimodal maps transform into families of graphs that turn out 
to be independent of map nonlinearity or other particulars. Here we provide an in depth description 
of the HV treatment of the Feigenbaum scenario, together with analytical derivations that relate to 
the degree distributions, mean distances, clustering coefficients, etc., associated to the bifurcation 
cascades and their accumulation points. We describe how the resultant families of graphs can 
be framed into a renormalization group scheme in which fixed-point graphs reveal their scaling 
properties. These fixed points are then re-derived from an entropy optimization process defined 
for the graph sets, confirming a suggested connection between renormalization group and entropy 
optimization. Finally, we provide analytical and numerical results for the graph entropy and show 

' that it emulates the Lyapunov exponent of the map independently of its sign. 

a \ 

'"O ! PACS numbers: 05.45.Tp, 05.45.Ac, 89.75.Hc 

o 

In recent years a new general framework to make time series analysis has been coined. This frame- 
work is based on the mapping of a time series into a network representation and the subsequent graph 
theoretical analysis of the network, offering the possibility of describing the structure of complex sig- 
nals and the associated dynamical systems from a new and complementary viewpoint, and with a 
full set of alternative measures. Here we focus on a specific type of mapping called the horizontal 
visibility algorithm, and via this approach we address the specific case of the period-doubling route to 

■ chaos. We extend our preliminary results on this topic [1], and provide a complete graph theoretical 
■^j- ' characterization of unimodal iterated maps undergoing period doubling route to chaos that, we show, 
t-H , evidence a universal character. Our approach allows us to visualize, classify and characterize periodic, 

' chaotic and onset of chaos dynamics in terms of their associated networks. 

<N ; 

O ' I. INTRODUCTION 

<N : 

Very recently a connection between nonlinear dynamical systems and complex networks has been accounted 
for by means of the horizontal visibility (HV) algorithm 0, as the latter transforms time series into graphs. The 
families of trajectories generated by nonlinear low-dimensional iterated maps conform a distinctive class of time series. 
Accordingly, they make up ideal candidates to test the capabilities of the HV algorithm for capturing meaningfully the 
information contained in them, and, if so, see how these manifest in the network central quantities. The possibility of 
observation of novel properties adds to the motivation to carry on these studies. We have chosen to inspect first the 
well-known one-dimensional (therefore dissipative) unimodal maps and their common period-doubling route to chaos, 
when periodic attractors transform into aperiodic attractors, the bifurcation cascade or Feigenbaum scenario [HQ. 
This route to chaos appears an infinite number of times amongst the family of attractors generated by unimodal maps 
within the windows of periodic attractors that interrupt sections of chaotic attractors. In the opposite direction, a 
route out of chaos accompanies each period-doubling cascade by a chaotic band-splitting cascade, and their shared 
bifurcation accumulation points form transitions between order and chaos that possess universal properties 043- 
Low-dimensional dynamics benefits from added interest as systems with many degrees of freedom relevant to various 
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problems in physics and elsewhere are known to undergo a drastic simplification and display this type of dynamics 
0- 

There is a growing number of methods designed to transform series into networks, involving concepts such as 
recurrence in phase space or Markov processes [llj to cite a few, and our approach forms part of this enterprise 

Q . Once a time series is converted into a network the interest lies in the observation of the characteristic properties 
of dynamical systems in a different environment. And to achieve this it is necessary to use the characteristic tools 
of network analysis fl2liil| . The family of visibility algorithms has been successful in obtaining information relevant 
to the description of fractal behavior [TtJ or to the distinction between random and chaotic series Q . Here we detail 
the Feigenbaum scenario as seen through the HV formalism by providing a complete description of its associated 
set of graphs. These graphs represent the time evolution of all trajectories that take place within the attractors of 
unimodal maps. The outline of the presentation is the following: We start in Section II by recalling the construction 
of an HV graph from a time series and deduce general expressions for the mean degree and distance when the 
series is periodic. We advance a visual illustration of the graphs and their location in the Feigenbaum diagram. 
In Section III we focus on the period-doubling cascade and derive a simple closed-form expression for the degree 
distribution of periodic attractor graphs and their accumulation point. We obtain from the latter the mean degree, 
the variance and the clustering coefficient. In Section IV we center on the reverse bifurcation cascade of chaotic-band 
attractors and derive the expression for the degree distribution and the mean degree. As the number of bands 
increases there is a growing similarity with the same quantities for the period-doubling cascade, since, as shown, the 
contribution from chaotic motion is confined only to the shrinking top band. The properties of graphs stemming 
from both chaotic-band attractors and windows of periodicity are derived with help of the self-affine properties of 
the bifurcation cascades. In Section V we describe a renormalization group (RG) transformation, equivalent to the 
original functional composition RG transformation, but specially designed for the Feigenbaum graphs, that leads 
to a set of fixed point graphs that further explain, and give unity to, the two previous sections. In Section VI we 
turn attention to the entropy associated to the degree distributions and find that under optimization we recover the 
RG fixed points. Finally we compare the behavior of this entropy as we move along the bifurcation cascades and 
notice that this quantity follows closely the variation of the map's Lyapunov exponent, pointing out to a property 
reminiscent of the Pesin equality but suitable for both periodic and chaotic graphs. In Section VII we summarize our 
results. A brief preliminary account of the contents of this paper is given in 



II. FEIGENBAUM GRAPHS 

The HV graph Q associated with a given time series {xi}i=i,...,Ar of N real data is constructed as follows: First, a 
node i is assigned to each datum Xi, and then two nodes i and j are connected if the corresponding data fulfill the 
criterion Xi, Xj > x n for all n such that i < n < j. Let us now focus on the Logistic map [H defined by the quadratic 
difference equation x t +i — f(x t ) = fJ,x t (l — x t ) where x t £ [0,1] and the control parameter fx £ [0,4]. According 
to the HV algorithm, a time series generated by the Logistic map for a specific value of its control parameter /x 
(after an initial transient of approach to the attractor) is converted into a Feigenbaum graph (see figured]). Notice 
that this is a well-defined subclass of HV graphs where consecutive nodes of degree k = 2, that is, consecutive data 
with the same value, do not appear, what is actually the case for series extracted from maps (besides the trivial 
case of a constant series). Also, as proven in (l8| . an HV graph is, by construction, a planar graph, that is, it has 
a diagram representation in which any pair of links intersect only at their endpoints. Moreover, an HV graph is 
also outerplanar: each node contacts the infinite face, where a face is a bounded region of a planar graph and the 
infinite face is its outer region. In what follows we take advantage of these facts and outline some generic properties 
of Feigenbaum graphs. 

Mean degree k. Consider a periodic orbit of period T. Without lack of generality, we represent the orbit as 
the infinite time scries {..., xo, xi, Xt, %x, x%, ...}, where xq = xt corresponds to the largest value of the series. By 
construction, the associated Feigenbaum graph consists of a concatenation of identical motifs of T+ 1 nodes associated 
to the subseries {xq, x\, xt}- Suppose that the motif is a graph with V links, and let xi be the smallest datum of 
the subseries which, by construction, will have degree k — 2 (and since no data repetitions are allowed in the motif, 
Xi will always be well defined). Now remove this node and its two links from the motif. The resulting motif will have 
V — 2 links and T nodes. Iterate this operation T — 1 times (see figure [5] for a graphical illustration of this process in 
a particular case with T = 5). The resulting graph will have only two nodes, associated with xq and xt, connected 
by a single link, and the total number of deleted links will be 2(T — 1). The mean degree k of the graph corresponds 
to the mean degree of the motif made of T nodes (the nodes associated with xq and xt only introduce half of their 
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FIG. 1: Feigenbaum graphs from the Logistic map Xt+i = f(xt) = fiXt(l — Xt). The main figure portrays the family of attractors 
of the Logistic map and indicates a transition from periodic to chaotic behavior at fioo = 3.569946... through period-doubling 
bifurcations. For n > //oo the figure shows merging of chaotic-band attractors where aperiodic behavior appears interrupted by 
windows that, when entered from their left-hand side, display periodic motion of period T = m ■ 2° with m > 1 (for /i < /too, 
m — 1) that subsequently develops into m period-doubling cascades with new accumulation points /ioo(?n). Each accumulation 
point Hooirn) is in turn the limit of a chaotic-band reverse bifurcation cascade with m initial chaotic bands, reminiscent of 
the self-afnne structure of the entire diagram. All unimodal maps exhibit a period-doubling route to chaos with universal 
asymptotic scaling ratios between successive bifurcations that depend only on the order of the nonlinearity of the map [IsJ, the 
Logistic map belongs to the quadratic case. Adjoining the main figure, we show time series and their associated Feigenbaum 
graphs according to the HV mapping criterion for several values of /i where the map evidences both regular and chaotic behavior 
(see the text). Inset: Numerical values of the mean normalized distance d as a function of mean degree k of the Feigenbaum 
graphs for 3 < fi < 4 (associated to time series of 1500 data after a transient and a step Sfi — 0.05), in good agreement with 
the theoretical linear relation (see the text). 



degree in the motif, what is equivalent to an effective reduction of one node). Hence, 

r , . Redoes 2(2(T - 1) + 1) -, s ( 1 \ 

The above result holds for every periodic or aperiodic (T — > oo) series, independent of the deterministic process 
that generates them, as the only constraint in its derivation is that data within a period are not repeated. It 
therefore includes all graphs generated by unimodal maps irrespective of their degree of nonlinearity. Observe that 
the maximum mean degree (achieved for aperiodic series) is k(oo) — 4, in agreement with previous theory (see fl8jV 

Normalized mean distance d. On the other hand, the normalized mean distance d of the graph is defined as 
d = DJN, where D is the mean distance (the average over all pairs of nodes of the smallest path that connects each 
pair) and N the number of nodes. For graphs associated with periodic orbits d depends on T (as this is the maximal 
amount of nodes that can be jumped through a link), and straightforwardly gives d(T) = ^ for N —> oo. Therefore, 
for HV graphs d and k are linearly related by 



J(fc) (4_fc). 



(2) 
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FIG. 2: Graphical illustration of the constructive proof of the expression for the mean degree k(T) through consideration of 
a motif extracted from a periodic series of period T = 5. Observe that the second data is the one with the lowest value. By 
deleting it, the graph loses 2 links. This process is iterated three more times. A total amount of 2 ■ (5 — 1) = 8 links has been 
deleted, independently of the inner structure of the T = 5 motif. 



This latter analytical relation is checked numerically in the inset of figure [T] The limiting solution k — > 4, d — > 
holds for all aperiodic, chaotic or random series. In addition to the numerical results shown in the inset of figure [T] for 
the specific case of the Logistic map, we have also examined the accuracy of the latter relation for several unimodal 
maps, giving perfect agreement in every case (data not shown). 



III. PERIOD-DOUBLING ROUTE TO CHAOS: RESULTS 
A. Order of visits of stable branches and chaotic bands 



A deep-seated feature of the period-doubling cascade is that the order in which the positions of a periodic attractor 
are visited is universal [20j . That is, the visiting order of the positions {xi} (where the subindex i denotes the 
iteration time) of a periodic attractor along the period-doubling route to chaos is the same for all unimodal maps 
20]. This ordering turns out to be a decisive property in the derivation of the structure of the Feigenbaum graphs 
(see figure [3] where we plot the graphs for a family of attractors of increasing period T = 2™, that is, for increasing 
values of /i < /Xoo)- Here we describe the rule that such ordering follows for orbits of period T = 2", and how this in 
turn induces the structure of the associated Feigenbaum graphs. This is illustrated graphically in figure [4] 

Consider the first period-doubling bifurcation that produces attractors with period T — 2 1 and for which repeated 
jumps are observed between two positions in time, x\, x%, X\, with x% > X2 (figure SJa). Without lack of generality 
label x\ as the largest data. This series is transformed into a Feigenbaum graph made up of a concatenation of a 
root motif of 3 nodes, where by construction the inner node is associated with datum X2- As the family of attractors 
reaches the next period-doubling bifurcation, each point Xi of the period-2 attractor splits into two new stable 
'offspring' positions: Xi and x^+t, and the visiting order is such that x± > a; 3 > x^ > X2 (figure S]b). This ordering is 



5 



n = 3 



n = 4 






o 

— 1 

CD 
£D 
W 

CO 

■p 



FIG. 3: Periodic Feigenbaum graphs for /i < /loo. The sequence of graphs associated to periodic attractors with increasing 
period T = 2 n undergoing a period-doubling cascade. The pattern that occurs for increasing values of the period is related 
to the universal ordering with which an orbit visits the points of the attractor. Observe that the hierarchical self-similarity of 
these graphs requires that the graph for n — 1 is a subgraph of that for n. 



reminiscent of the T — 2 1 orbit that was present before the bifurcation (namely, the orbit returns to a neighborhood 
of the point after a journey along the attractor). In particular, this means that the second largest value X3 (the 
offspring of xi) is visited only after a T = 2 1 journey, that is, in the middle of the T = 2 2 periodic orbit. Observe also 
that the bottom pair of offspring positions appears inverted (grey box in the figure). The corresponding Feigenbaum 
graph is a concatenation of the T = 2 1 motif (blue and red portions in the figure) linked by the largest node x\, which 
repeats after T = 2 2 iteration times. Blue and red portions in the graph are equivalent since the orbit follows the 
same pattern of visits across the stable branches (each portion according to a given offspring). The same procedure 
can be iterated for increasing period-doubling bifurcations (see figure HJc), leading to Feigenbaum graphs which are 
progressively self-similar and become so in the limit n — ¥ 00. Summing up, the period-doubling bifurcation of an orbit 
of period T — 2™ generates two identical copies of the T = 2™ root motif of the graph that are now concatenated by 
the node associated to datum X1+2™ and linked by the bounding nodes x%, and this in turn is the root motif of the 
T = 2™ +1 Feigenbaum graph. In the following section we will take advantage of this structure to analytically derive 
several topological properties of the Feigenbaum graphs along the period-doubling cascade. 



B. Topological properties of Feigenbaum graphs along the period-doubling cascade 

1, Degree distribution P(n,k) 

The above-described order of visits generates a hierarchy of self-similar Feigenbaum graphs along the period-doubling 
bifurcation cascade. The degree distribution of a graph is denned as a discrete probability distribution that expresses 
the probability of finding a node with degree k [13H15j . By construction, the degree distribution of a Feigenbaum 
graph for a series of period T — 2", n = 0, 1, 2, is 

P(n,k) = (5) 7 . fc = 2,4,6,...,2n, 
P(n,k)= , fc = 2(n+l), 

P(n,k)=Q, fc = 3,5,7, or fc>2(n+l). (3) 
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FIG. 4: Graphical illustration that explains how the order of visits to the stable branches of the map induces the structure of 
the Feigenbaum graphs all along the period-doubling bifurcation cascade (/i < jitoo)- 




FIG. 5: Dots: Semi-log plot of the degree distribution of a Feigenbaum graph associated with a time series of 10 data extracted 
from a Logistic map at the onset of chaos /ioo = 3.5699456.... The straight line corresponds to equation [4] in agreement with 
the numerical calculation (the deviation for large values of the degree are due to finite size effects). 

At the accumulation point = 3.5699456... the period diverges (n — > oo), and the degree distribution of the 
Feigenbaum graph at the onset of chaos becomes a (non truncated) exponential for even values of the degree, 

P(oo,A) = Q) i fc = 2,4,6,..., 

P(oc,fc)=0, fc = 3,5,7,... (4) 

In figure [5] we show the accuracy with which this analytical result is reproduced by a finite series of 10 6 data. 
Numerical and theoretical distributions are in good agreement. 
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2. Mean degree k(n) and normalized distance d(n) 
The mean degree k(n) is the first moment of the degree distribution 

2(»+l) , 1 v 

k(n) = J2 kP(n,k) = 4 (l- )' ( 5 ) 

fc=2 ^ ' 

which at the accumulation point yields fc(oo) — 4, indeed the maximal degree for an HV graph. 

Furthermore, by induction it can be shown that the mean distance D(n, N) within a Feigenbaum graph of N nodes 
after n period-doubling bifurcations is given approximately by 

D(n,N) « ^ |^±i + i [(6n - 1)2" + (-1)"]} . (6) 

Observe that for a fixed period T — 2™ the mean distance increases linearly with the size N of the graph. Normalization 
of this measure in the limit of infinite size leads to a well-defined mean distance per node d(n) with the simple form 

d(n) = ~, (7) 
such that the linear relation depicted in equation [5] is recovered. 



3. Clustering coefficient c(n, k) 

The local cluste ring coefficient is a topological measure that quantify how close a given node's neighbors are to 
being a clique fljMLa ]. The local clustering coefficient of a Feigenbaum graph that corresponds to an attractor of 
period T — 2" is given by 



k - 1 2 



c(n, k) = k = — , k — 2, 4, 6, 2n, 



V2/ 

c( "' 4) = ^ = |Kp ^ 2( " +1) ' (8) 

a result indicative of a so-called hierarchical structure [2l| . Since the degree distribution P(n, k) is known in closed 
form throughout the period-doubling cascade, use of the approximation c(n, 2(n + l))~l/(n+l) leads to 

, x fc/a 

P(n,fe)=f-J =P(n,2/c(n,fc)), fc = 2, 4, 6, 2n, 2(n + 1), 

P(n,c)= Qy^. (9) 
Consequently, the mean clustering coefficient c(n) is given by 

c-(„) . g c(„,«P(n, fc ) . g- (-) + ( „ +1)(2 „ +1) (5 

n , 

El . . 

m2 m + 2 n ~ 1 (n + l)(2n + 1) ' ^ ^ 

This last summation does not possess a solution in closed form except in the limit n — > oo which yields c(oo) = log 2 = 
0.693... Nevertheless, the series converges to log 2 extremely fast: c(0) = 0, c(l) = 0.666..., c(2) = 0.691666..., c(3) = 
0.693452..., which suggests that c(n) rapidly loses its dependence on the bifurcation order n and remains basically 
constant for all n. 



4- Higher moments of the degree distribution: variance o~ 2 (n) 



The moments of the degree distribution P(n, k) can be easily calculated by making use of the generating function 

M(t) = (e tk ) = l+U-^)t+ 



12- 



10 + 4n\ 2 /104 (100/3) + 20n + 4n 2 \ 3 , 4 



In particular, the variance is given by 



)t 3 + o (; 4 ) 



which at the accumulation point becomes cr 2 (oo) = 8. 



IV. REVERSE BIFURCATION CASCADE OF CHAOTIC BANDS: RESULTS 
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FIG. 6: Aperiodic Feigenbaum graphs for /i > pioo. A sequence of graphs associated with chaotic series after n chaotic-band 
reverse bifurcations, starting at fj. = 4 for n — 0, when the attractor extends along a single band and the degree distribution 
does not present any regularity (red links). For n > the phase space is partitioned in 2" disconnected chaotic bands and 
the n-th self-afhne image of = 4 is the n-th Misiurewicz point /i 2 „_i_ 2 „. In all cases, the orbit visits each chaotic band in 
the same order as in the periodic region fi < jitoo • This order of visits induces an ordered structure in the graphs (black links) 
analogous to that found for the period-doubling cascade. 



The Logistic map at ji = 4 is said to be fully chaotic as the Lyapunov exponent attains its maximum value of log 2 
there, and also the single-band attractor spans the unit interval. As [i is decreased smoothly from fi = 4 towards 
/loo, this first chaotic band [n — 0) suffers a contraction and a series of successive band splittings (see figure [1]) that 
are in several respects analogous to the period-doubling bifurcations at fi < fj,^. Since the band splittings occur 
from right to left, they are called reverse bifurcations [4TI0I [20|. After the first reverse bifurcation (n = 1), the first 
chaotic band splits into two disconnected chaotic bands, such that the iteration of an initial condition in this zone 
generates an orbit that alternates between these two bands. While the exact location of each position gives evidence 
of sensitive dependence on initial conditions the band alternation is fixed. Significantly, while in the chaotic zone 
orbits are aperiodic, for reasons of continuity they visit each of the 2™ chaotic bands in the same order as positions are 
visited in the attractors of period T = 2™ [20]. In figure [H] we have plotted the Feigenbaum graphs generated through 
chaotic time series at different values of jx that correspond to an increasing number of reverse bifurcations. 
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FIG. 7: Self-affinity in the chaotic region: the two disconnected chaotic bands at ^2-4 < ft < fJ>i-2 are rescaled copies of the 
first chaotic band < fJ> < 4. An orbit at /12-4 < fi < /ii-2 makes an alternating journey between both bands. 

The first of the reverse bifurcation points, called Misiurewicz points, is located at /.*i_2 = 3.6785.... As in the 
period-doubling region // < /loo, a cascade of chaotic band reverse bifurcations takes place in the complementary 
region /i > /i^, such that after n bifurcations the attractor consists of 2" non-overlapping chaotic bands. The 
location of the Misiurewicz points converges to //qo in the limit n — > 00, and their relative locations also obey the 
same asymptotic ratio 5 = 4.669... (Feigenbaum constant for quadratic maps) as n increases as that occurring for the 
period-doubling bifurcations at /i < /i,^ 0, [22j • 

Furthermore, the chaotic bands that are formed after each reverse bifurcation are self-affine copies of the previous 
ones [23| . This is illustrated in figure [7J where it can be observed that the two chaotic bands (blue and orange 
boxes) appearing after the first reverse bifurcation (at /ii-2) are rescaled copies of the first chaotic band (red box). 
Accordingly, each value of [i within the first chaotic band has n self-affine images after n reverse bifurcations. For 
instance, the first Misiurewicz point /ii-2 = 3.6785... is a self-affine image of fi = 4 after the first reverse bifurcation, 
the second Misiurewicz point /i2-4 = 3.5925... is also a self-affine image of \i = 4 after the second reverse bifurcation, 
and so on. Consider now a chaotic time series extracted at /i G [^2-4 — ^1-2] in which the iterate position values 
alternate between both chaotic bands. Now separate this series into two time-ordered subseries, each of which contain 
only data belonging to either the top or the bottom chaotic bands. These subseries are indeed rescaled copies of a 
series extracted from the first chaotic band at the corresponding preimage value of fi. Notice that the HV algorithm 
is invariant under affine transformations in the series. Accordingly, the Feigenbaum graphs of these subseries are the 
same as the graph obtained for the self-affine value of fi belonging to the first chaotic band. 

A. Self-affine properties of chaotic bands: mean degree and degree distribution 

As mentioned, the interband motion that an orbit experiences within an attractor composed of 2™ chaotic bands 
follows the same order of visits as that we described for a periodic orbit of period 2™ . And as in the case of the periodic 
region, this property possibilitates the derivation of the main characteristics of the Feigenbaum graphs in the chaotic 
region. Consider the Feigenbaum graph of a chaotic series generated by the Logistic map after n reverse bifurcations. 
By construction, the structure of this graph is in many respects a mirror image of the Feigenbaum graph in the 
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periodic region after n period-doubling bifurcations. The only differences originate from the data located within the 
top chaotic band. This feature is due to the fact that chaotic bands do not overlap and that the HV algorithm filters 
out the precise locations of data in favor of unspecified relative positions. Denote by P M (n, k) its degree distribution, 
where, in order to not overburden the notation, we refer to fi indistinctly of the self-affine image of fi we study. By 
construction, and up to k = 2n, the degree distribution P M (n, k) is equivalent to its periodic region counterpart, i.e. 

fl\ k/2 

PMin,k)=l-\ , fc = 2,4,6,...,2n. (11) 

Now, the data belonging to the top chaotic band transmits a specific contribution to the associated Feigenbaum 
graph which in general cannot be determined analytically, and as we shall see shows up in P M (n, k) for k > In. 
This contribution is denoted P^ op (n, k), where, by definition, P^ op (0, k) — P M (0,fc), the degree distribution of the 
Feigenbaum graph at the first chaotic band n = 0, fii-2 < fJ- < 4. While the precise shape of the ^-dependent 
distribution Pff p (n,k) is unknown, the self-affine structure of the band-splitting cascade described previously makes 
it possible to relate it to the first chaotic band, i.e. 

P*°P(n,fc)=P M (0,fe). 

This important property, that, as we shall comment below, corresponds to a crossover phenomenon in the RG flows, 
breaks down the structure of the Feigenbaum graphs within the chaotic region into two contributions: (i) an interband 
contribution stemming from the order of visits amongst chaotic bands that is equivalent to that for periodic attractors 
(with associated black links in figure [5]) , and (ii) a contribution from the first chaotic band (with associated red links 
in figure [6]). For normalization reasons, this second contribution can be coarse-grained as a cumulative distribution, 
leading to 

P M (n,fc) - (^j ' , fc = 2,4,6,...,2n, 

P>,fc>2(n + 1))= . (12) 

Interestingly, as n increases the contributions from the interband motion of the chaotic attractor these become 
more and more dominant and the contribution associated with the first chaotic band is progressively smeared out, 
disappearing at the accumulation point /Xoo (see figure [5]). 

Based on the previous arguments on self- affinity, the contribution from the first chaotic band can be determined 
quantitatively as follows. For concreteness, let us focus on a chaotic attractor located after the first reverse bifurcation. 
The attractor is divided into two non-overlapping chaotic bands, and every orbit in it alternates between the 'top' 
and 'bottom' bands. While the precise location of the visits within each band is unknown, by construction the nodes 
associated to the data located in the 'bottom' band have degree k = 2, and this occurs for half of the iteration times, 
so 

W>2) = \- 

Also, for parity reasons (even number of chaotic bands) one has 

P„(l,3) = 0, 

independently of fi. For k > 3, the contribution comes necessarily from the top band, and normalization yields 

P^(l,fc) = ^P* op (l.fc-2), fc>4. 

Finally, self-affinity implies 

P^(l,k)=P fl (0,k). 

Repeating this same procedure, we find that after n reverse bifurcations the degree distribution of a Feigenbaum 
graph fulfills the expressions 



fl\ k ' 2 

P Al (n,fc)=f-J , k= 2,4,6,..., 2n, 



(13) 
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P fl (n,k) = 0, k = 3,5,7, ...,2n+ 1, (14) 

independently of fj,, and 

P»(n,k)= P;°P(n,fc-2n)= P M (0,fc-2n), fc>2(n + l), (15) 

that depends on As expected, the last expression only depends on the structure of the first chaotic band, decreases 
for increasing values of n and disappears in the limit n — > oo. 

The above expressions for P /i (n, k) can be used to derive the mean degree kp,{n) of a Feigenbaum graph in the chaotic 
regime: 



2/ ^ pv ' V 2™ / 2 n 

fc=2 v 7 x 7 fe=2(n+l) 

This last expression relates the mean degree after n reverse bifurcations to that in the first chaotic band. Interestingly, 
the only solution with constant mean degree in the chaotic regime is k^O) = kn(n) = 4, in full agreement with the 
general theory. 



B. Periodic windows: self-affine copies of the Feingenbaum diagram 

The Feigenbaum diagram shows a rich self-affine structure: periodic windows of initial period m that undergo 
successive period-doubling bifurcations with new accumulation points /ioo(w) appear interwoven with chaotic 
attractors at \x > /i^. These period-doubling cascades taking place beyond fj,oo are self-affine copies of the 
fundamental period-doubling route to chaos (/i < /loo)- For instance, the window that initiates with period m = 3 
(/i w 3.84) generates three period-doubling cascades, each one being a properly rescaled copy of the fundamental 
one (see figure [5] for a graphical illustration) . In order to take into account such phenomenology in the labeling 
of Feigenbaum graphs, we make use of the following notation: G(m, n) is the Feigenbaum graph associated with 
a periodic attractor of period T — m • 2 n , that is, the graph belonging to a periodic window of initial period m, 
after n period-doubling bifurcations. Observe that the process of chaotic band reverse bifurcations that take place 
from fj, = 4 towards /i^ is again repeated in a self-affine manner, indeed, each accumulation point /ioo(m) is the 
limiting value of a chaotic-band reverse bifurcation cascade. In figure |8] this aspect is also illustrated. Accordingly, 
we may also extend the notation of Feigenbaum graphs belonging to chaotic regions, such that G^fm^n) is associ- 
ated with a chaotic attractor composed by m-2™ bands (that is, after n reverse bifurcations of m initial chaotic bands). 

Given a periodic window with initial period m within the first chaotic band (fJ>i-2 < M < 4), we note that after n 
reverse bifurcations, the chaotic-band attractors contain a self-affine image of that window, where there is an initial 
periodic attractor of period m ■ 2™ developing into m • 2™ period-doubling cascades inside the window. In this case, 
equation [16] reads 

M „, m . 2 -) = 4 (i-i) + M^, 

that together with equation [T] yields 

kJn,m- 2 n ) = 4(1- -— L— V 

M ' 1 \ 2(m2™)7' 

as expected. 



V. RENORMALIZATION GROUP APPROACH 

A. RG transformation: definition, flows and fixed points 

In order to recast previous findings in the context of the renormalization group, let us define an RG operation 1Z 
on a graph as the coarse-graining of every couple of adjacent nodes where one of them has degree k = 2 into a block 
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2 9 3 3.1 3.! 3 3 3 4 3.5 3.B 3.7 3.8 3.9 4.0 3.330 3 835 3 340 3 845 3.850 3.8S6 

U is 

FIG. 8: Zoom of the Feigenbaum diagram close to the period m = 3 window. Starting with a period 3 orbit, each one of the 
stable branches develops into a period-doubling bifurcation cascade with a new accumulation point fj,oo(3), beyond which the 
attractor becomes chaotic, interwoven with periodic windows: each part of the diagram is indeed a rescaled copy of the full 
Feigenbaum tree. The locations of several Feigenbaum graphs (with the notation defined in the text) within the period three 
window are depicted. 



node that inherits the links of the previous two nodes (see figure [S]a). This is a real-space RG transformation on the 
Feigenb aum graph [24], dissimilar from recently suggested box-covering complex network renormalization schemes 



271 ] . As a matter of fact, this scheme turns out to be equivalent for /i < /loo to the construction of an HV graph 



from the composed map instead of the original /, in correspondence to the original Feigenbaum renormalization 
procedure [1, [Hj]. We first note that JZ{G(l,n)} — (7(1, n — 1), thus, an iteration of this process yields an RG flow 
that converges to the (1st) trivial fixed point TZ^{G(1, n)} = (2(1,0) = Go = TZ{Gq}- This is the stable fixed point 
of the RG flow V/i < /ioo. We note that there is only one relevant variable in our RG scheme, represented by the 
reduced control parameter A/i = /ioo — hence, to identify a nontrivial fixed point we set A/i = or equivalently 
ri — > oo, where the structure of the Feigenbaum graph turns to be completely self-similar under 1Z. Therefore we 
conclude that G(l, oo) = Goo is the nontrivial fixed point of the RG flow, lZ{Goo} — Goo- In connection with this, let 
Pt(k) be the degree distribution of a generic Feigenbaum graph Gt in the period-doubling cascade after t iterations of 
1Z, and point out that the RG operation, 1Z{Gt} = Gt+i, implies a recurrence relation (1 — P t (2))P t +i{k) — P t (k + 2), 
whose fixed point coincides with the degree distribution found in equation 2) This confirms that the nontrivial fixed 
point of the flow is indeed Goo ■ 



Next, under the same RG transformation, the self-affine structure of the family of chaotic attractors yields 
TZ{Gfj,(l,n)} — G^(l,n — 1), generating a RG flow that converges to the Feigenbaum graph associated to the 1st 
chaotic band, TZS n ^{G^(l, n)} = G^(1,0). Repeated application of 1Z breaks temporal correlations in the series, 
and the RG flow leads to a 2nd trivial fixed point TZ^°°^{G^(l, 0)} = G ranc i = 7?.{G ran d}, where G ran d is the 
HV graph generated by a purely uncorrelated random process. This graph has a universal degree distribution 
P(k) — (l/3)(2/3) fc_2 , independent of the random process underlying probability density (see @, Q). 

Finally, let us consider the RG flow inside a given periodic window of initial period m. As the renormalization 
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FIG. 9: Renormalization process and network RG flow structure, (a) Illustration of the renormalization process 1Z: a node 
with degree k — 2 is coarse-grained with one of its neighbors (indistinctively) into a block node that inherits the links of both 
nodes. This process coarse-grains every node with degree k — 2 present at each renormalization step, (b) Example of an 
iterated renormalization process in a sample Feigenbaum graph at a periodic window with initial period m — 9 after n — 2 
period-doubling bifurcations (an orbit of period T — m-2" — 36). (c) RG flow diagram, where m identifies the periodic window 
that is initiated with period m and n designates the order of the bifurcation, n — n + 1 for period-doubling bifurcations and 
n — — (n + 1) for reverse bifurcations. A/i(m) = /loo(fi) — /i denotes the reduced control parameter of the map, and (J.oo(m) 
is the location of the accumulation point of the bifurcation cascades within that window. Feigenbaum graphs associated with 
periodic series (A/i(m) > 0, n > 0) converge to G(1,0) = Go under the RG, whereas those associated with aperiodic ones 
(A/i(m) < 0, n < 0) converge to G ra nd- The accumulation point fioo = M°o(l) corresponds to the unstable (nontrivial) fixed 
point G(l, oo ) = Gx> of the RG flow, which is nonetheless approached through the critical manifold of graphs G(m, oo) at 
the accumulation points /Xoo (jn) ■ In summary, the nontrivial fixed point of the RG flow is only reached via the family of the 
accumulation points, otherwise the flow converges to trivial fixed points for periodic or chaotic regions. 

process addresses nodes with degree k = 2, the initial applications of 1Z only change the core structure of the graph 
associated with the specific value m (see figure [Sib for an illustrative example). The RG flow will therefore converge 
to the 1st trivial fixed point via the initial path TZ^{G(m, n)} — G(l,n), with p < m, whereas it converges to the 
2nd trivial fixed point for G M (m,n) via 7^^ p - ) {G At (m, n)} = G^{l,n). In the limit of n — s> oo the RG flow proceeds 
towards the nontrivial fixed point via the path TZ^' {G(m, oo)} = (7(1, oo). Incidentally, extending the definition of 
the reduced control parameter to A/i(m) = /ioo(w) — the family of accumulation points is found at A/i(m) = 0. A 
complete schematic representation of the RG flows can be seen in figure Oc. 

Interestingly, and at odds with standard RG applications to (asymptotically) scale-invariant systems, we find that 
invariance at A/i = is associated in this instance to an exponential (rather than power-law) function of the 
observables, concretely, that for the degree distribution. The reason is straightforward: 1Z is not a conformal 
transformation (i.e. a scale operation) as in the typical RG, but rather, a translation procedure. The associated 
invariant functions are therefore non homogeneous (with the property g(ax) = bg(x)), but exponential (with the 
property g(x + a) = cg(x)). 



B. Crossover phenomenon 

When the RG transformation for a Feigenbaum graph is applied repeatedly, it generates flows terminating at two 
different trivial fixed points (7(1, 0) = Go and G ran d or at a nontrivial fixed point G(l, oo) = Goo. Go is a chain graph 
where every node has two links, G ran d is a graph associated with a purely random uncorrelated process, whereas Goo is 
a self-similar graph that represents the onset of chaos. The RG properties within the periodic windows are incorporated 
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into a general RG flow diagram (see figure[9lc and ngure[TO]for an alternative representation). Here we add a comment 
on the standard presence of crossover phenomena in the RG applications, in our case for large n (or /i ~ /i^) for 
both fj, < fioo and fj, > fioo. In both cases the graphs (7(1, n — j) and G M (1, n — j) with j n closely resemble the 
self-similar G^ (obtained only when fi = /i^) for a range of values of the number j of repeated applications of the 
transformation TZ until a clear departure takes place towards Go or G ran( j when j becomes comparable to n. Hence, 
for instance the graph {G ^,(1 , n)} will only show its true chaotic nature (and therefore converge to G ran d) once 
j and n are of the same order. In other words, this happens once its degree distribution becomes dominated by the 
contribution of P^° p (n, k) (alternatively, once the graph core -related to the chaotic band structure and the order of 
visits to chaotic bands- is removed by the iteration of the renormalization process) . 




FIG. 10: Illustrative cartoon incorporating the RG flow of Feigenbaum graphs along the entire Feigenbaum diagram: aperiodic 
(chaotic or random) series generate graphs whose RG flow converges to the trivial fixed point G ra nd, whereas periodic series 
(both in the region /i < jitoo and inside periodic windows) generate graphs whose RG flow converges to the trivial fixed point 
G(0, 1). The nontrivial fixed point of the RG flow G(oo, 1) is only reached through the critical manifold of graphs at the 
accumulation points fioo(m). 



VI. NETWORK ENTROPY 

A. Entropy optimization and RG fixed points 

Recently [28], it has been point out that there exists a connection between the extremal properties of entropy 
expressions and the renormalization group approach, when applied to systems with scaling symmetry. Namely, that 
the fixed points of RG flows can be obtained through a process of entropy optimization, providing the RG approach 
a variational flavor. In this section we investigate on this connection, and derive, via optimization of an entropic 
functional for the Feigenbaum graphs, all the RG flow directions and fixed points directly from the information 
contained in the degree distribution. We depart by applying Jayne's principle of Maximum Entropy [29] with two 
different restrictions and prove that the degree distribution P(k) that maximizes its entropy h is in each case related 
to the fixed points of the RG flow. We employ the standard optimization technique that involves Lagrange multipliers. 

Let us begin by defining a graph entropy as the Shannon entropy of the graph degree distribution 



oo 

h=-Y / P(k)logP(k), 

k=2 



(17) 
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and assume that the degree distribution P(k) has a well-defined mean k (this is true in general for HV graphs and in 
particular for Feigenbaum graphs according to equation [1]). The, consider the Lagrangian 

oo / oo \ / oo \ 

£ = -^P(fe)logP(A : )-(Ao-l) (£p(fc)-l) -A x ^feP(fc)-fc , 

fc=2 \k=2 / \fc=2 / 

for which the extremum condition reads 

dC 



logP(fc) - A - Ai/c = 0, 



dP{k) 

and has the general solution 

P{k) = e - x "- Xlk . 

The Lagrange multipliers Ao and Ai can be calculated from their associated constraints. First, the normalization of 
the probability density, 

oo 



k=2 

implies the following relation between Aq and Ai 



oc 



e A i - r 

fc=2 

and differentiation of this last expression with respect to Ai yields 

e- Xl - 2 



fc=2 V ; 

Second, the assumption that the mean degree is a well-defined quantity (true for HV graphs) yields 



S^ke- X °- Xlk = k = 
^ 1 - e~ Xl 

k=2 



Combining the above results we find 



Ai = log 



and 

A = lo. 

Hence, the degree distribution that maximizes h is 

P(k) = 



k-1 
k-2 



k-1 



k-1 (k-2 xk 



(k - 2) 2 \k - 1 

which is an increasing function of k. The maximal entropy is therefore found for the maximal mean degree k = 4, 
with an associated degree distribution 

3 /2\ fc 1 /2 X fe 2 



P(fc) 4 V 3y 3 V 3 

Remarkably, we conclude that the HV graph with maximal entropy is that associated with a purely uncorrelated 
random process, and coincides with G lan d, a trivial fixed point of the RG flow. 
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Observe now that, by construction, the Feigenbaum graphs along the period-doubling route to chaos (// < /loo) do 
not have odd values for the degree. Let us assume this additional constraint in the former entropy optimization pro- 
cedure. The derivation proceeds along similar steps, although summations now run only over even terms. Concretely, 
we have 



oo 1 

E -Ai2fe 1 
e 2A! _ l ■ 



k=l 



which after differentiation over Ai gives 



e 



2Ai _ 2 



k=l 



(gSAi _ 1) 



2 



and 



We obtain for the Lagrange multipliers 



and 



^2fce- A °- Al2fe = k 



2e 



2Ai 



fc=l 



=2Ai _ l ' 



Ai = - log 



2 ° V fc - 2 



A = log 



fc-2 



2 

The degree distribution that maximizes the graph entropy turns now to be 



P(k) = 



2 (k-2^ k/2 



fc - 2 



As before, entropy is an increasing function of fc, attaining its larger value for the upper bound value k = 4, which 
reduces to P{k) = (l/2) fe ^ 2 , fc = 2,4,6,..., equation|3J We conclude that the maximum entropy of the entire family 
of Feigenbaum graphs, if we require that odd values for the degree are not allowed, is achieved at the accumulation 
point, that is, the nontrivial fixed point of the RG flow. 

Finally, the network entropy is trivially minimized for a degree distribution P(2) = 1, that is, at the trivial fixed point 
Go • These results indicate that the fixed point structure of a RG flow can be obtained from an entropy optimization 
process, confirming the aforementioned connection. 



B. Network entropy and Pesin theorem 

1. Periodic attractors 

By making use of the expression we have for the degree distribution P[n, fc) in the region \i < (j.^ we obtain for 
the graph entropy h(n), after the n-th period-doubling bifurcation, the following result 

2(n+l) 

h(n) = - P(n,k) log P(n, k) 

k=2 

= -E^ lo g(^)-^iog(^)=^(fc-|r) = iog4(i-^). (18) 

We observe that the entropy increases with n and, interestingly, depends linearly on the mean degree fc. This 
linear dependence between h and fc is related to the fact that, generally speaking, the entropy and the mean of a 
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probability distribution are proportional for exponentially distributed functions, a property that holds exactly in 
the accumulation point (eqHJ) and approximately in the periodic region (eq[3]), where there is a finite cut off that 
ends the exponential law. Interestingly, HV graphs associated with chaotic series have also been found to have an 
asymptotically exponential degree distribution [3(. Finally, note that in the limit n — > oo (accumulation point) the 
entropy converges to a finite value h(oo) — log 4. 



2. Chaotic attractors 

For Feigenbaum graphs (7^(1, n) (in the chaotic region), in general h cannot be derived exactly since the precise 
shape of P(k) is unknown (albeit the asymptotic shape is also exponential Q). However, arguments of self-affinity 
similar to those used for describing the degree distribution of Feigenbaum graphs can be employed in order to find 
some regularity properties of the entropy h^in). Concretely, the entropy after n chaotic band reverse bifurcations 
can be expressed as a function of n and of the entropy in the first chaotic band h^(0). Using the expression for the 
degree distribution, a little algebra yields 



h^n) = log 4 + 



log 4 - 




FIG. 11: Horizontal visibility network entropy h and Lyapunov exponent A for the Logistic map. We plot the numerical values 
of h and A for 3.5 < fi < 4 (the numerical step is Sfi = 5 • 1CP 4 and in each case the processed time series have a size of 2 12 
data). The inset reproduces the same data but with a rescaled entropy h — log(4). The surprisingly good match between both 
quantities is reminiscent of the Pesin identity (see text). Unexpectedly, the Lyapunov exponent within the periodic windows 
(A < inside the chaotic region) is also well captured by h. 



In figure [T] we observe how the chaotic-band reverse bifurcation process takes place in the chaotic region from 
right to left, and therefore leads in this case to a decrease of entropy with an asymptotic value of log 4 for n — > oo at 
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the accumulation point. These results suggest that the graph entropy behaves qualitatively as the map's Lyapunov 
exponent A, with the peculiarity of having a shift of log 4, as confirmed numerically in figure [TT1 This unexpected 
qualitative agreement is reasonable in the chaotic region in view of the Pesin theorem that relates the positive 
Lyapunov exponents of a map with its Kolmogorov-Sinai entropy (akin to a topological entropy) that for unimodal 
maps reads hxs = A, VA > 0, since h can be understood as a proxy for hxs- Unexpectedly, this qualitative 
agreement seems also valid in the periodic windows (A < 0), since the graph entropy is positive and approximately 
varies with the value of the associated (negative) Lyapunov exponent even though hxs = 0, hinting at a Pesin- like 
relation valid also out of chaos which deserves further investigation. The agreement between both quantities lead us 
to conclude that the Feigenbaum graphs capture not only the period-doubling route to chaos in a universal way, but 
also inherits the main feature of chaos, i.e. sensitivity to initial conditions. 



VII. SUMMARY 

In summary, we have described how the horizontal visibility algorithm transforms nonlinear dynamics into networks, 
in particular how the entire Feigenbaum scenario manifests in network language. The most important quality of the 
HV algorithm when applied to this category of time series is that it leads to definite and unique results, to well-defined 
families of graphs and to analytic expressions for their main property, the degree distribution. In obtaining these 
results two basic ingredients played a decisive role; these are the predetermined order in which positions or bands 
are visited in the time series and the self-affinity of periodic or chaotic-band attractors. As we have seen, the graphs 
generated by the algorithm store in the connectivity pattern amongst nodes, i.e. in their topological structure, the 
dynamical nature of the nonlinear map. A significant feature of the Feigenbaum graphs is their evident universality as 
the structure of all graphs and that of all families they form is independent of the specific form of the unimodal map 
including the degree of nonlinearity of its extremum. We have also shown that the families of networks and degree 
distributions obtained from periodic and chaotic attractor bifurcation cascades have scale-invariant limiting forms. 
And that the latter occupy the dominant positions of RG fixed points and extrema of the entropy associated with the 
degree distribution. The capability of the HV algorithm to expose new information is indicated by the property of 
the network entropy to emulate the Lyapunov exponent for both periodic and chaotic attractors. Applications of this 
approach to other low-dimensional nonlinear circumstances, such as the dynamical complexity associated with the 
quasiperiodic route to chaos and the vanishing of Lyapunov exponents are under study, and an obvious open question 
is how well would these precise analytical results work for experimental data. 
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